
ELSEVIER 


Available online at www.sciencedirect.conn - 

ScienceDirect PHYSICS LETTERS A 

Physics Letters A 372 (2008) 3423-3428 

www.elsevier.com/locate/pla 


Multiple branches of ordered states of polymer ensembles 
with the Onsager excluded volume potential 

Hongyun Wang® *, Hong Zhou'’ 

^ Department of Applied Mathematics and Statistics, University of California, Santa Cruz, CA 95064, USA 
^ Department of Applied Mathematics, Naval Postgraduate School, Monterey, CA 93943, USA 

Received 4 December 2007; received in revised form 22 January 2008; accepted 24 January 2008 
Available online 16 February 2008 
Communicated by C.R. Doering 


Abstract 

We study the branches of equilibrium states of rigid polymer rods with the Onsager excluded volume potential in two-dimensional space. Since 
the probability density and the potential are related by the Boltzmann relation at equilibrium, we represent an equilibrium state using the Fourier 
coefficients of the Onsager potential. We derive a non-linear system for the Fourier coefficients of the equilibrium state. We describe a procedure 
for solving the non-linear system. The procedure yields multiple branches of ordered states. This suggests that the phase diagram of rigid polymer 
rods with the Onsager potential has a more complex structure than that with the Maier-Saupe potential. A study of free energy indicates that 
the first branch of ordered states is stable while the subsequent branches are unstable. However, the instability of the subsequent branches does 
not mean they are not interesting. Each of these unstable branches, under certain external potential, can be made metastable, and thus may be 
observed. 

© 2008 Elsevier B.V. All rights reserved. 
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The isotropic-nematic phase transition in rigid rod-like poly¬ 
mers is a classical topic [9]. In 1949, Onsager first examined 
the isotropic-nematic phase transition theoretically by model¬ 
ing the steric excluded volume interactions [1]. Onsager’s study 
was based on a virial expansion, which yields a mean field po¬ 
tential, now bearing his name (the Onsager potential). Using 
this mean field potential and choosing an appropriate trial func¬ 
tion for the orientation distribution function, Onsager was able 
to argue that when the concentration is high enough, there is 
a transition from a uniform isotropic state to an orientationally 
ordered prolate nematic state. Many theories have been formu¬ 
lated after Onsager’s pioneering work. An excellent review on 
the theoretical advances of liquid crystalline polymers is given 
in [2]. In this study, we adopt the Doi-Hess model [3,4], which 
was first developed by Doi and Edwards to describe spatially 
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homogeneous flows of rodlike liquid crystal polymers and has 
been used in many studies [5-7]. In the Doi-Hess model, the 
polymer ensemble is represented by an orientational probabil¬ 
ity density in a meso-scale much larger than individual polymer 
rods but much smaller than the macroscopic flow. Each poly¬ 
mer rod undergoes Brownian diffusion and is affected by the 
interaction with other polymer rods within the meso-scale. The 
inter-molecular interaction is usually modeled using the On¬ 
sager potential or the Maier-Saupe potential, which is an ap¬ 
proximation to the Onsager potential. The main mathematical 
difference between the Onsager potential and the Maier-Saupe 
potential is that the Onsager potential depends on the whole 
probability density function while the Maier-Saupe potential 
depends only on the second moment. As a result, the Doi-Hess 
model with the Onsager potential is mathematically much more 
challenging than with the Maier-Saupe potential. 

The Doi-Hess model with the Onsager potential has been 
studied numerically and the Onsager theory has been extended 
to other systems. Eor example, Larson [8] applied the spherical 
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harmonic expansions to solve the Doi-Hess equation with the 
Onsager potential for the 3D time-dependent orientation distri¬ 
bution function in the presence of shearing flow. Lasher [10] 
extended Onsager’s work to the nematic ordering of hard rods 
with a direct application of the scaled particle approach. Vroege 
and Lekkerkerker [11] provided a comprehensive overview on 
the theory and experiments in lyotropic colloidal and poly¬ 
mer liquid crystals. In particular, they generalized Onsager’s 
virial theory to poly disperse solutions and soft interactions. 
Chrzanowska [12] developed a simple Onsager theory type den¬ 
sity functional theory (DFT) of a two-dimensional system of 
hard needles. 

Recently, the phase diagram of rigid rod polymers with 
the Maier-Saupe potential, previously obtained in various nu¬ 
merical studies (for example, see [5] and [13] and references 
therein), has been studied analytically [14-16]. Thus, it is 
timely and worthwhile to revisit analytically the phase diagram 
of rigid rod polymers with the Onsager potential. In studying a 
phase diagram, it is important to get a full picture of the phase 
diagram, including all stable and unstable branches. Besides 
mathematical completeness, there are two more reasons for in¬ 
cluding all branches when studying a phase diagram: (a) the 
full structure of the phase diagram helps us understand the dy¬ 
namics and the stabilities; and (b) under a certain external field, 
unstable branches may be made meta-stable and thus may be 
observed. 

In this Letter, we consider the case of rigid rod polymers in 
the two-dimensional space with the Onsager intermolecular in¬ 
teraction potential. In the two-dimensional space, orientation is 
represented by the polar angle 0 and the orientational probabil¬ 
ity density p(0,t) is governed by the Smoluchowski equation 


/ ^Onsager -h —^ 


( 1 ) 


where D is the rotational diffusion coefficient of each polymer 
rod, ks the Boltzmann constant and T the absolute temperature. 
For simplicity, we use U{0) to denote the Onsager potential 
normalized by which has the form 
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= b j \s,in0\p{0 + 0)d0. (2) 

0 

Here b is proportional to the normalized polymer concentration 
and is inversely proportional to the temperature. We start with 
the Fourier expansion of probability density p{6) 
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The Fourier expansion of p(0 -h 0) as a function of 0 has the 
form 

^ oo 

p (0 -h 0) =-h / (ak cos kS -h bk sin kS) cos kS 
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Substituting Fourier expansion (6) into Onsager potential (2), 
and using identities (4) and (5), we obtain the Fourier expansion 
of 1/(0) 
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The equilibrium probability density is related to the Onsager 
potential by the Boltzmann distribution: 


p(0) = 2exp(-Ui0)), 

27T 

Z = j exp{-U{0))d0. (8) 

0 

The Boltzmann relation allows us to use the Fourier coefficients 
of 1/(0) to represent the equilibrium state. For that purpose, we 
introduce 
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In terms of {c 2 j,d 2 j}, U(S) and p(0) have the expressions: 


We first express the Onsager potential in terms of the Fourier 
coefficients {ak, bk} of p(0). For that goal, we first derive two 
integral identities. 


U(S) = 


40 

In 


oo 

y^(c 27 cos2j0 -h d 2 j sin2j0), 
;=i 


( 10 ) 









H. Wang, H. Zhou/Physics Utters A 372 (2008) 3423-3428 


3425 


p{9) = —exp{-U(e)) 


^ exp ^ ^(C2 j cos 2j9 + d 2 j sin 2j9)^ 


( 11 ) 


where Z is the partition function given by 

Z = y* expf cos 2j0 + d 2 j sin2j0)jd0. 

Eq. (11) shows that the equilibrium probability density p(0) is 
periodic with period tt . Consequently, in the Fourier expansion 
of p( 0 ), we have a 2 j-\-i = ^ 2 ;+i = 0 - 

Next we show that U(0) can be made an even function by 
shifting the coordinate system. Using (9) to relate p(0) given in 
(3) to the derivatives of U(0) given in (10), we get 

I ^ 
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Combining this result with the Boltzmann relation ( 8 ) leads to a 
non-linear differential equation for the Onsager potential U{0)\ 

[/"(0) + [/(0)-yexp(-[/(0))=O. (13) 

Let ^0 be the location where U(0) attains its maximum in 
[0, 27r]. We shift the coordinate system to make = 0- Thus, 
the initial conditions for differential equation (13) at 0 = 0 are 
U (0) = Uq and U'iQi) = 0. Both the differential equation and 
the initial conditions are invariant under the change of variable 
6 ^^^ = —^oid- Therefore, we conclude that U{0) = U{— 6 ), and 
consequently d 2 j = 0. It follows from (9) that b 2 j = 0. In terms 
of {c 2 j}, U( 0 ) and p( 0 ) have the expressions: 
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(9) is one set of equations relating Fourier coefficients {a 2 j } and 
{c 2 ;}. Calculating Fourier coefficients {a 2 j} directly from p(0) 
expressed in (15) yields a second set of equations relating {a 2 j} 
and {C 2 j]: 
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Combining (9) and (16), we arrive at a non-linear system for 
{C 2 k}’ 

cos 2k0 exp(^^^ C 2 j cos2 jO)dO 
Jo"" exp(J]~ 1 C 2 j cos2j9)d9 
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Below we describe a procedure for solving non-linear sys¬ 
tem (17). Here, we use the Fourier coefficients {c 2 ;} of U(0) 
to represent an equilibrium state. The equilibrium probability 
density p(0) is related to U{0) by the Boltzmann distribution 
and is expressed in terms of {c 2 ;} in (15). The advantage of 
using Fourier coefficients {c 2 ;} and the exponential form (15) 
is that the positivity of equilibrium probability density p( 0 ) is 
automatically enforced without putting any constraints on the 
Fourier coefficients {c 2 ;}. 

We solve non-linear system (17) for ordered states. Here an 
ordered state is one in which not all the coefficients of {c 2 ;} are 
zero. We first find ordered states with C 2 0. In non-linear sys¬ 
tem (17), Z? is a physical parameter that can be tuned by varying 
the polymer concentration or the temperature. From the point 
view of physics, it is desirable to express (c 2 , C 4 , c^,...), the 
solution of nonlinear system (17), as a function of b. Mathemat¬ 
ically, however, it is more convenient to treat C 2 as an indepen¬ 
dent variable and view Z? as a function of C 2 . We adopt this math¬ 
ematical formulation. Specifically, in non-linear system (17), 
we let r = C 2 , we treat r as a parameter, and we solve for b 
and (c 4 , C 6 , C 8 ,.. .)• The iterative procedure described below 
provides a way of constructing solution b and (C 4 , C 6 , eg,...) 
for a given value of C 2 = r 7 ^ 0. For conciseness of presenta¬ 
tion, we introduce shorthand notations. Let C = (C 4 , C6, eg,.. .)• 
Remember C 2 is fixed at the given value of r. We introduce 
functions: 


cos 2 k 0 exp(r cos 26 -h YIJL 2 ^ 


exp(r cos 20 -h ^ 2 ; cos 2 j 0 )d 0 

1,2,3,.... (18) 

With these notations, non-linear system (17) becomes 
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Cg^...). Remember that, in our mathematical formulation, 
C 2 = r 7 ^ 0 is given while b is part of the solution to be de¬ 
termined. So we use Eq. (19) to update b and then use Eq. (20) 
to update the rest of the Fourier coefficients (C 4 , C 6 , cg,...). 


Update b: : 
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Our iterative procedure ^ ^(^+1) is 
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This is how we calculate the first branch of ordered state. The 
convergence of iterative method (21) is verified in numerical 
implementations. A rigorous mathematical proof for the con¬ 
vergence, however, is still open. 

Next we calculate the second branch of ordered states by 
finding solutions of (17) with C2 = 0 but C4 ^ 0. Specifically, 
in non-linear system (17), we keep C2 = 0, we keep C4 7^ 0 
fixed and treat it as a parameter, and we solve for b and 
(c6, eg, CIO,.. .)• One way to keep C2 = 0 is to assume that so¬ 
lutions on the second branch have 2x2 = 4 fold symmetry, 
that is, they are periodic with period 27r/4. If g(^) is periodic 
with period 27r/4, then we always have cos 20g{0) dO = 0. 
Our numerical calculations indicate that invoking symmetry is 
the only way to keep C2 = 0, although a rigorous mathemati¬ 
cal proof is still open. With the 4 fold symmetry, we only need 
to consider the Fourier coefficients {c 4 j, j ^ 1}. Recall that in 
our mathematical formulation, we keep the leading non-zero 
Fourier coefficient fixed and treat it as a parameter. For the sec¬ 
ond branch, we keep C4 7^ 0 fixed and treat it as a parameter. 
Then, the resulting non-linear system for b and { 04 j, j > 1} 
is solved using an iterative procedure similar to (21) described 
above. The subsequent branches of ordered states are calculated 
in a similar way. For example, the third branch is calculated by 
invoking 2x3 = 6 fold symmetry. 

To plot the phase diagram, we represent an ordered state 

using the magnitude of the Onsager potential, • 

Fig. 1(a) shows the first three branches of ordered states of 
rigid rod polymers with the Onsager potential. The first branch 
(labeled as Si{b)) intersects with the isotropic state at Z? ~ 
4.7124. The second branch (labeled as S 2 {b)) intersects with 
the isotropic state at b ~ 23.5619. The third branch (labeled 
as S^Qy)) intersects with the isotropic state at Z? ~ 54.9779. 
The procedure of calculating ordered states can be contin¬ 
ued beyond the first three branches shown. Numerical results 
suggest that there are infinitely many branches of ordered 
states where the kth branch is characterized by having 2 k- 
fold symmetry with C 2 k = 0 but C4y^ 7^ 0. For this reason, in 
the space of {c2;, j = 1 ,..., 00}, which is of infinite dimen¬ 
sions, two different branches of ordered states will never inter¬ 
sect each other. Interestingly, in the two-dimensional space of 

as shown in Fig. 1(a), it is also true that 
two different branches of ordered states will not intersect each 
other. 

In comparison, in the case of polymers with the Maier- 
Saupe intermolecular potential, there is only one branch of or¬ 
dered states and it intersects with the isotropic state at Z? = 4. 
Fig. 1(b) compares the only branch of ordered states for the 
Maier-Saupe polymer (dashed line) with the first branch of or¬ 
dered states for the Onsager polymers (solid line). 



Fig. 1. (a) The first three branches of ordered state for the Onsager polymer, 
(b) The comparison of the only branch of ordered states for the Maier-Saupe 
polymer (dashed line) with the first branch of ordered states for the Onsager 
polymers (solid line). 


In the numerical calculation above, we used {c2;, j = 
1,..., A} with N = 4096. The infinite sum •jYlJLi(.^ 2 j)^ is 
approximated by the corresponding partial sum. To show that 
this approach is at least self-consistent and that N = 4096 is 
large enough, we examine how fast |c2; | converges to zero as j 
increases. Fig. 2(a) shows the plot of |c2; | vs index j for the so¬ 
lution on the first branch at Z? = 80. It is clear that as j increases, 
|c2; I decays exponentially to zero, and the limit of machine pre¬ 
cision is reached at about j ~ 500. 

To determine the stability of the branches of ordered states 
shown in Fig. 1(a), we examine the free energies of these equi¬ 
librium states. The Onsager free energy functional of probabil¬ 
ity density p(0) at normalized polymer concentration b is given 
by 


2n 

G{[p\b)^ j (^p{e)inp{e) + ^u{e)p{e)ye. ( 22 ) 
0 

Fig. 2(b) shows the free energies of the isotropic branch 
(dashed line) and the first three branches of ordered states (solid 
lines) for the Onsager polymers. More detailed studies of free 
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(a) 


lytically so that all conclusions are established on solid mathe¬ 
matical ground. Once the 2D case of the Onsager polymers is 
analytically solved, we will move on to revisit the 3D case an¬ 
alytically. For the 3D case of the Onsager polymers, we shall 
start with a rigorous mathematical derivation of the axisym- 
metry. The axisymmetry of equilibria of the 3D Maier-Saupe 
polymers has been mathematically established only recently in 
three studies [14-16]. All these three studies rely heavily on the 
specific form of the Maier-Saupe potential. For analytical stud¬ 
ies of the 3D Onsager polymers, new mathematical tools are 
needed. 

The multiple branches of ordered states for the Onsager 
polymers stem from the fact that the Fourier expansion of func¬ 
tion I sin 01, the kernel in the Onsager potential, has infinite 
number of terms. 



0 20 40 , 60 80 100 

b 

(b) 

Fig. 2. (a) Plot of \c 2 j\ vs index j. (b) Free energies of the isotropic branch 
(dashed line) and the first three branches of ordered states (solid lines) for the 
Onsager polymers. 


energies reveal that the isotropic branch is stable before the first 
branch of ordered states appears (i.e., for b < 4.7124). Once 
the first branch of ordered states appears (i.e., for b > 4.7124), 
the isotropic branch becomes unstable. The first branch of or¬ 
dered states is stable while the subsequent branches ordered 
states are unstable. The instability of the subsequent branches 
of ordered states does not mean these branches are not inter¬ 
esting. For example, the second branch of ordered states can 
be made metastable under an external potential of the form 
d^Externai(^) = A COS 4^, and thus may be experimentally ob¬ 
served under such a condition. 

In [17], the equilibrium bifurcation map of 3D polymers 
with the Onsager interaction potential was studied. A linear sta¬ 
bility analysis of the base isotropic state was carried out. It was 
found that high order bifurcations are linearly unstable. Our re¬ 
sults on 2D polymers with the Onsager interaction potential 
are consistent with the 3D results in [17]. In a future study, 
we will (a) tackle the mathematical conjectures suggested by 
numerical calculations; and (b) carry out asymptotic analysis 
to identify the exact bifurcation locations and find asymptotic 
expansions of branches near bifurcation locations. One of our 
goals is to first re-derive all conclusions in the 2D case ana- 


I sin^l 


2 

TT 


00 


E 


4 

7t{4k^ — 1 ) 


cos2k0. 


(23) 


The solutions on the kth branch have 2k fold symmetry. Basi¬ 
cally, in every meso-scale neighborhood (since we are studying 
spatially homogeneous solutions) the ensemble of polymer rods 
is divided into k sub-groups of equal size, and each sub-group 
is aligned with a director. The k directors are uniformly spaced 
on the unit circle, resulting in a 2k fold symmetry. This 2k fold 
symmetry of the kth branch corresponds to the terms cos 2k0 , 
cos 4k0 and beyond in the Onsager potential (23). Indeed, for an 
ensemble of polymer rods with the Onsager interaction poten¬ 
tial, if k sub-groups and k directors are sustained by an external 
force, then a newly added polymer rod will align with any of 
these k directors individually with equal probability. Again, 
a newly added polymer rod is able to see these k directors in¬ 
dividually because of the terms cos2k0, cos4k0 and beyond 
in the Onsager potential (23). In the Maier-Saupe potential, 
however, the kernel is replaced by — cos^ 0 = const — ^ cos 20 , 
which has only the cos 20 term in its Fourier expansion (ex¬ 
cluding the constant term). As a result, the 2k fold symmetry 
is suppressed for k > 1. For an ensemble of polymer rods with 
the Maier-Saupe interaction potential, even if k sub-groups and 
k directors are sustained by an external force, a newly added 
polymer rod will not be able to see the k directors individu¬ 
ally for k > 1. That is why the polymers with the Maier-Saupe 
potential have only one branch of ordered states. 
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